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A nonlinear model for the vibration suppression of a smart 
composite elastic plate using graphical representation 
involving fuzzy control is presented. The plate follows the 
von Karman and Kirchhoff plate bending theories and the 
oscillations are caused by external transversal loading forces, 
which are applied directly on it. Two different control forces, 
one continuous and one located at discrete points, are 
considered. The mechanical model is spatially discretized by 
using the time spectral Galerkin and collocation methods. 
The aim is to suppress vibrations through a_ simulation 
process within a modern graphical computing environment. 
Here we use MATLAB/SIMULINK, while other similar 
packages can be used as well. The nonlinear controller is 
designed, based on an application of a Mamdani-type fuzzy 
inference system. A computational algorithm, proposed and 
tested here is not only effective but robust as well. 
Furthermore, all elements of the study can be replaced or 
extended, due to the flexibility of the used SIMULINK 
environment. 


1. Introduction 


Classical mathematical theories of control work well on linear systems. However, their 
effectiveness depends on many restrictions. On the other hand, nonlinear controllers, based on 
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fuzzy logic with built-in smart computational methods can give a good description of a behavior 
of nonlinear structures and provide nonlinear feedback. The details of the comparison between 
the classical control and fuzzy logic approach are discussed, among others, in the paper of 
Tairidis et al. [1]. 


Active vibration control of smart elastic plates has been considered by many authors. An active 
vibration control incorporating active piezoelectric actuator and self-learning control for a 
flexible plate structure is described by Tavakolpour et al. [2]. In the paper of Fisco and Adeli [3] 
hybrid control strategies, in particular, active and semi-active vibration suppression are 
discussed. A brief survey on industrial applications of fuzzy control in different fields of 
engineering is given in the works of Precup and Hellendoorn [4]. The control systems are 
classified into three groups, control systems with Mamdani fuzzy controllers, control systems 
with Takagi-Sugeno fuzzy controllers, and adaptive and predictive control systems. A review of 
active structural control is done by Korkmaz in [5]. 


Here a vibration suppression of a smart plate is performed by the algorithm, involving 
constructed model-diagrams of SIMULINK with the use of Mamdani FIS. In the paper [6] a 
Newmark technique is applied for solving the system of linear ordinary differential equations 
(ODEs), obtained after the spatial discretization of the model. In this paper, we use a graphical 
representation for solving linear and nonlinear ODEs. The numerical simulation is implemented 
with the use of model-diagrams of SIMULINK. That allows easily and quickly suppressing 
vibrations of the plate and improves the presentation of the results. We can also include a 
Sugeno-type fuzzy inference system in the complete model-diagram for the linear model, 1.¢e., the 
new approach provides flexibility of choosing a controller. Besides, we have tested a collocation 
method, which provides an opportunity to locate control forces at discrete points of the plate. 
The proposed algorithm allows effective suppression of the Fourier coefficients, which display 
characteristics of displacement, velocity and Airy's stress function (in other words, general 
displacement, velocity and Airy's function) of the plate. After computing the coefficients, we can 
easily calculate the displacement and velocity at each point of the plate. 


Active fuzzy control is a suitable tool for the systematic development of nonlinear control 
strategies and can be fine-tuned if no experience exists in the behavior of the considered body 
(structure) or if one designs more complicated control schemes, (e.g. [7—11]). Similar approaches 
for active vibration control of a simply supported rectangular plate using fuzzy logic rules are 
discussed in the paper of Shirazi et al. [11]. The results have been compared with the classical 
proportional integral derivative (PID) controller. 


The present paper is organized as follows. Section 2 focuses on a formulation of the nonlinear 
dynamic mechanical model and set up initial and boundary conditions. In Section 3 the spatially 
discretized model is presented. In both cases (Galerkin' projections and collocation approach) the 
result of discretization is a system of second-order nonlinear ordinary differential equations of 
motion. In Section 4 an algorithm for vibration suppression is introduced. The algorithm 
involves, constructed model-diagrams with composed “S-Function’” of SIMULINK. A numerical 
example using Galerkin's projections and the developed algorithm is presented. Section 5 is 
devoted to a state space representation of the linearized spatially discretized problem. We 
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reformulate algorithm mentioned above for the linear model. Instead of the *‘S-Function" a 
““State-Space" block mode' is activated. The implementation of the collocation method for the 
linear case is illustrated in a numerical example. The main results are reported in Section 6. 


2. Formulation of the problem 


Governing equations. A nonlinear mathematical model describing the bending vibrations of an 
isotropic, homogeneous elastic plate in the presence of active control, and allowance for the 
rotational inertia of the plate elements and viscous damping is written as 


L,(w,W) = phwy — p™ Awe + hew, + DA*w — h[w, p] = [Q] + [Z], (1) 
L.(w,p) = Ap -< [w,w]=0, (t.x,y)€Q. (2) 


Here the following notations are used: 
[w, y] = 04, W022 + 01,W022.W = 2012W012W (Monge-Ampére’s form). 


o wis the deflection (displacement) of the plate. 

o wW(t,x,y) is the Airy stress potential describing internal stresses, which appear due to the 
deformation of the plate (e.g., [12—-14]). 

o 2 = (0,T| X G, where T is the final time and G = (0,1,) x (0, 1,) is the shape of the plate (1, 

and l, are the lengths of the sides of the plate). 

p is the density of the material. 

h is the thickness of the plate. 

c is the viscous damping coefficient. 

D is the flexural rigidity of the plate. 

[Q| are the external transversal loading forces. 

[Z| are the control forces. 


00000 0 


The equations Eq.1, Eq.2 are extensions of the nonlinear von Karman plate model for large 
deflections ( [13,14], etc.) on a control case. The plate is subjected to external transversal 
disturbances and generalized control forces, produced, for example, by electromechanical 
coupling effects. 


For the Galerkin’s projections, it is required that w, y € C?(0,T; W**(G)), (W” is the Sobolev 
space). For the collocation method, it is required that w, p € C?(0,T;C?(G)) U C(0,T; C*(G)). 


Regarding the loading forces, we consider [Q] as a function of (t,x, y) or as a discrete function 
only of t, defined at some collocation points of the plate. We assume the same for the control 
forces [Z]. When [Z] is supposed to be a function of (t, x, y), it is expanded into double Fourier's 
series and we use the time spectral Galerkin method for spatial discretization of L,, Lz. The 
suppression of vibrations is done through the control of the Fourier coefficients. Here we also 
introduce a new approach, by considering a time-discrete type of control force [Z], located at 
some discrete-collocation points of the plate, which may be different from the points where 
external forces are applied. 
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Initial and boundary conditions. For the displacement and velocity at the initial time, it is 
assumed: 
w(0,x,y) = u(x,y), w,(0,x, y) = v(x,y) in G, 
where the functions u,v belong to the quadratically integrable function space L2(G). 


Regarding the boundary conditions, we consider simply supported, partially and clamped plates 
(see [15,16]). 


3. Spatial discretization of the problem Eq.1, Eq.2 


3.1. Time spectral expansions 


An approximate analytical solution of Eq.1, Eq.2 in the form of partial sums of double Fourier's 
series with the time-dependent coefficients ([15]) reads 


N N 
Wye.) = > wi Ooyy), Welter y=) wih © wy), 3) 
i,j=1 ijat 
ul ne 
WEN= > wWOgyay),  (t>0), (4) 
ij=1 


where the global basis functions w,; are chosen to match the boundary conditions [15]. For the 
initial conditions we have 


N 
Wy (0.2.9) = uv Gey) = >. wil (0) @yGe¥), 5) 
ow pe 
We (0,29) = wu(xy) =) wi (0) oy C6) (6) 
ij=1 


where @;; and oj; are the bases. It is assumed jj; = 0); = @;;. In the case when [Z] is a function 
of (t,x, y) we expand it into double Fourier series with the basic functions Wj, i.¢., 


N 
IZ) = 2y(txy) = >. 2h) wy(x,y). (7) 
ij=l 


According to the two algorithms, introduced first for the linear case in [6], and extended for the 
nonlinear case in [16] the vibration suppression can be performed by means of the control 


function Z or the Fourier coefficients zi (t). In the present paper, the second case is studied. 


A. Muradova et al./ Journal of Soft Computing in Civil Engineering 2-1 (2018) 01-17 5 


3.2. Galerkin's projections 


Let the external loading force [Q] is a quadratically integrable function Q, i.e., Q € L2(G) and 
[Z] is defined as Eq.7. Introducing the inner product in Lz - space and applying Galerkin's 
projections to Eq.1, Eq.2 yields 


(Ly (Wy, Py ), mn) = (Q, Omn) + (Zn, Omn); 
(L>(Wy, Py), mn) =0, mn =1,2,...,N. 
Hence 
Mwy (t) + CWy(t) + Kywy(t) = Ary(wy(t), Wy (t)) + Pan (t) + Fzy (0), (8) 
K,Wy(t) = Azan (wy (t), wy (t)), (9) 


where Wy(t), Wy (t) are the vectors with the components which are the Fourier coefficients for 
the displacement Eq.3, and the Airy stress function Eq.4, respectively. We call them 
characteristics of the displacement and the Airy stress function, respectively. Further, M = 
ph(H + (h?/12)B) is the mass matrix (B is the approximation of the Laplacian), C is the 
viscous damping matrix, K, is the stiffness matrix. The matrices K,, Kz are the result of 
approximation of the biharmonic operator and Aj, Az. are nonlinear approximations of the 
Monge-Ampére forms. Finally, P and F are the matrices and qy, Zy are the vectors, obtained 
after the approximations of the external exciting pressure and control forces, respectively. The 
operators H, B, P and F take different forms depending on the boundary conditions. The 
description of the above defined operators are given in [15]. Substituting Eq.9 into Eq.8 we 
obtain 


MWy + CWy + KyWy = Aay (Wy, K2*A2,v (Wy, Wy) + Pqw + Faw. (10) 


For the components from Eq.10 we have 


(MWy mn + (CWy) mn + (KiWy) mn = (4. (ww, K>*A2 (Wy, wy))) 


mn 
+(Pqu)mn + (FZu)mn mn = 1,2,..,N. (11) 


For the initial conditions Eq.5, Eq.6 after applying Galerkin's projections we obtain 


(Uwy (0), = (Un (%Y), Omn)s 
(Vwy(0)) = (Uy(%,Y),@mnn), mn=1,2,...,N. 
where u, v are described in [15]. 


If the inputs for the fuzzy controller are wy'"(t) and wy'"(t) then the output is the control 
Zmn(t) for these coefficients. 
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3.3. Collocation approach 


Let us now consider the collocation points (xz, v1), {(z, V1), O< x, <4, O< y, <b, k,l= 
1,2,...,N} on the spatial domain G for the system Eq.1, Eq. 2. Suppose that [Q] is the function 
of time defined at some or all collocation points. Similary let the control forces [Z] be put at 
some collocation points, which may coincide with the points of application of the loading 
pressure. We find the solution of Eq.1, Eq.2 in the form Eq.3, Eq.4 when Wy and Wy satisfy the 
equations Eq.1, Eq.2 at the collocation points, i.e. 


L,(Wy, Pw =xpy=y1) = [Ql w=xpy=¥7) a [ZJleaxpv=v1)) (12) 
L2 (Wy, Uw lexpyan) = 9 t= 1,2,..,N. (13) 


Supposing that control forces are located at every or some collocation points, from Eq.12 and 
Eq.13, we obtain a system of nonlinear ordinary differential equations of motion with respect to 
wy” and wi”, similar with the previous one Eq.10 


— — — — —-1— 
Mwy + CWy + Ki Wy = Ain (ws K, RCW, Wy) + q + Z y (14) 


and similar to Eq.11 


ma ae wy a = 1— —_— _ 
(Mwy), + (Cww),, + (Kiwy),, = (Ran (ws K, Ran(¥n 9) + Gi + Zl > 
kl 
k,l=1,2,...,N, (15) 


where M is the mass matrix, C is the damping matrix and K, is the stiffness matrix. The 
elements of the matrices are determined through computations of the partial derivatives of Wy 
and Wy (see Eq.3, Eq.4) with respect to the spatial variables x,y at the collocation points. The 
elements of the matrices M , C , K, and K, can be easily calculated. They are the values of the 
basic functions w;; (xx, y,). Analogously, the nonlinear parts are defined. 


Furthermore, in Eq.14 q is a vector with components q,,, the values of the time-discrete forces at 
some collocation points (xx, y,), k,l = N,,N,+1,...,N2 and Z is the vector with components 
Zz, the values of the control forces at some collocation points (x,,y,), k,l = M,,M,+ 
1,...,M 2. . The collocation points where the external loading forces are applied are called 
loading points and the points where we put/locate the control are called control points. In case 
N, = M,< N and N,= M,< N the loading points coincide with the control points. 
Obviously, we can also deal with free collocation points, where the external and control forces 
are absent. At these points, the values of the external and control forces are supposed to be zero. 


The advantage of the collocation method over the Galerkin's projections is that we do not need to 
take the inner products, and we can consider external loading disturbances at the discrete points 
and locate the control forces at the places of (all or some) collocation points in the proper way. 
Inversely, the collocation points can be considered at the best positions for the control, 1.e., the 
collocation points are chosen in order to provide optimal suppressions of vibrations of the plate. 
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A disadvantage of the collocation method is that one may lose accuracy of the approximate 
solution between the collocation points, and the obtained mass, damping and stiffness matrices 
are not so sparse as in case of Galerkin's projections. 


For the initial conditions from Eq.5, Eq.6 we have 


N N 
>, wh wy.) = www, —-D. wely (Oeoij.y) = vv te WO, 


i,j=1 j=l 


3.4. State-space representation 


Let us now suppose v,; = WwW, V2 =W (here and below the index N is omitted for the 
convenience). Then from Eq.10, we obtain 


V1 = V> 
Mv, + Cv, + Kv; = Ain (v., Kz*A2 y(Vz, v,)) + Pqy + FZy 
or 
Vi = V2, 
V> = M1 |-cv, = KV; + Pqy + FZy + Ain (vi, Kz*A2 y (V1, v:))}. (16) 


Analogously, for the model Eq.14 
Vi = V2, 
V> = M 


—, — — ——-1i1— 
|-c. = Kv + q + Z + Ain (v. K, Az (V4, v2] . (17) 


4. Algorithm of suppression of the Fourier coefficients and their derivatives 


The use of model-diagrams of SIMULINK allows us quickly and effectively calculate the 
suppressed Fourier coefficients in the expansions Eq.3, Eq.4. After that one can easily calculate 
the suppressed vibrations at each point of the plate. 


The decision of the controller is passed on the transformation (mapping) of input and output 
variables into membership functions along with a set of ‘‘rules". A total number of 15 rules are 
used here (Table 1). All rules have weight equal 1 and *‘AND" type logical operator is used. 


Table 1 
The fuzzy inference rules (e.g., if the displacement is “‘FarUp" and the velocity is ‘‘UP" then the control 


force is ‘‘Max") 
Vel\Disp. FarUp CloseUp Equil | CloseDn FarDn 
Up Max Med+ Low+ Null Low- 
Null Med+ Low+ Null Low- Med- 
Down Low+ Null Low- Med- Min 


8 A. Muradova et al./ Journal of Soft Computing in Civil Engineering 2-1 (2018) 01-17 


The Mamdani FIS has two inputs and one output. The inputs for the FIS are the displacement 
and velocity, and the output is the control force. Triangular and trapezoidal shape membership 
functions are chosen both for the inputs, Figures 1, 2 and for the output, Figure 3. 


input 1 MF 


se, = Equillbrum Close, 


Degree of menbestip 


-3 -2 = 0 
Displacement x10 


Fig. 1. Membership functions for the input (displacement) of the fuzzy inference system. 


nput 2 MF 


08 


o 
a 


Degree of menbeshp 
o 
ty 


02 


Velocity x10" 


Fig. 2. Membership functions for the input (velocity) of the fuzzy inference system. 
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Output MF 


Degree of menbeship 


Control Force 


Fig. 3. Membership functions for the output (control forces) of the fuzzy inference system. 


The implication method is set to minimum (min), and aggregation method is set to maximum 
(max). The defuzzified output value is created by using the Mean of the Maximum 
defuzzification method. 


A trial-and-error adjustment is used for tuning the coefficients of the membership functions. A 
fully automatic fine-tuning using genetic and particle swarm optimization ([17,18]) and also 
training the data with ANFIS ([19,20]) is possible and very effective for particular cases. 


Below a step by step algorithm for suppression of the Fourier coefficients in a SIMULINK- 
MatLab environment with the use of a Mamdani type, fuzzy inference system (FIS) is presented. 
For the nonlinear system two diagrams of SIMULINK, an initial one (Figure 4) and a complete 
model diagram (Figure 5) are employed in the simulation. The diagrams include **S-Function" 
option which allows the insertion of MatLab code into a SIMULINK block. In Figures 4, 5 “‘S- 
Function" is the MatLab code Sfun_NonLin. The initial module is intended for computations of 
the coefficients without control, and the complete module makes suppression of the coefficients. 
The suppression is performed by a nonlinear Mamdani fuzzy controller ([6,16]). 


In Figure 5 the suppression of the coefficients w’*, wz* in Eq.3 are provided on each time step 
of the simulation by the model-diagram. The control is put in the first equation in Eq.11 and 
Eq.15 for the first component, i.e. z‘t # 0 and the others are zero. The control function Eq.7 is 
easily estimated at every point (x,y). 
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Algorithm for vibration suppression 
Step 1. Set up input data of the initial diagram (Figure 4) in the block “Sine Wave" by defining 


an amplitude, frequency phase, etc. of exciting loading harmonic function Q(t) = p sin(wt + 
b) + @. Here p is the amplitude, w is the frequency, b is the bias and @ is the phase. 


Sfun_NonLi 


: (t) 
Sine Wave S-Function Plots 


Fig. 4. The initial model-diagram for the nonlinear system 


Step 2. Set up a number of basic functions N (formulas Eq.3, Eq.4). Correspondingly, determine 
the sizes of input and output data in the initialization part of ‘‘S-Function", initialize the 
initial conditions for the system Eq.15 and the array of sample time. 

Step 3. Run the initial module (without control) for obtaining a numerical solution of Eq.16 or 
Eq.17, i.e., in order to obtain the necessary data for the training of the controller. The 


output is the discrete values w°/(t;,) and, wi °(t,.) which are written as xi. 

Step 4. Generate Mamdani-type fuzzy inference system (FIS) using the results obtained in Step 
3, get the rules and then load the composed FIS on Workspace using Fuzzy Logic 
Toolbox: fuzzy. 

Step 5. Run the complete module, Figure 5 (*Selector" block corresponds to Example 1), with 


the use of the obtained FIS in order to get suppressed w‘0/0(t,,) and wo °(t,). (We can 


use this FIS for computations of the other coefficients, w/(t;), we (th), i # ip, j # jo). 
The complete model-diagram incorporates two subsystems (Subsystems 1, 2"). 
“Subsystem 1" is intended for getting results without control and “Subsystem 2" 
involves the control. The subsystems have the same structure with the initial model- 
diagram, presented in Figure 4. 
Step 6. Switch on computing another coefficient, namely change the outputs in “S-Function”, 
and go to Step 3. Otherwise, go to Step 1 or Step 2 or terminate the simulation. 
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Fig. 5. The complete model-diagram. 


In order to calculate the suppressed displacement W(t,,x,y) and suppressed velocity 
W, (ty, x,y) at any point (x,y) of the plate at discrete time t = t, , we use the values of the 
coefficients, obtained by the above given algorithm and substitute them into the analytical 
expressions Eq.3. The coefficients Wy (t) in Eq.4 can be computed through the equations Eq.9 or 
Eq.13. We can compute the coefficients as many as necessary in order to get a more accurate 
solution in sense of suppression of the displacement W(t, x, y) and velocity W;(t, x,y) at each 
point of the plate. 


Below we consider an example of vibration suppressions of a simply supported plate by the 
described above algorithm with the use of the time spectral (Galerkin) method, Section 3. 


Example 1. 

Let the external force be an exciting harmonic load, Q(t) = sin10z t uniformly distributed on 
the simply supported plate. Then Pynn = (Q, mn) = (P Gu) mn = 8/ih, /m*mn. The amplitude 
in the option for the *‘Sine Wave" block in Figures, 4, 5 is a vector P with 18 components. It is 
defined as P = (0,p11,0, P12, 0,713, 0, P21, 0, P22, 0, P23, 0, P31, 0, P32, 0, P33)". The control 
forces are added only to the second component p,,, i.e., to the right hand side of the second 
equation Eq.16. Thus, we have three ‘‘Selector" blocks: “1 of 18", “2 of 18" and “3-18 of 18". 
The control forces are added in the *‘Selector" block **2 of 18", 1.e., to the second component of 
the vector P. 
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For the physical parameters of the plate we have l, =l, =1, EF=2, v=0.3, h=0.5, p= 
1000, c=0, D =1. Furthermore, the number of basic functions N = 3 and the final time 
T = 30. The ODE solver chooses a time step with respect to the given error, initial and minimal 
time steps. The first coefficients in Eq.3 are calculated. The results are presented in Figures 6 and 
7 with the use of the Scope of SIMULINK. In Figure 6 we display the characteristics of the 
displacement and velocity, i.e., the first Fourier coefficients w1* and w;' before and after the 
control. In Figure 7 the external initial and control forces are shown, respectively. The x-axis 
(horizontal) denotes t - the time of simulation in seconds while the vertical axes denote the first 
Fourier coefficients in the expansions for the displacement and velocity, Eq.3. 


Displacement (coef) before (green --) and after (red —) flurry control 


x 10 
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_ 


Displacement (coef) 
i - ° 


0 5 10 16 20 26 wn 
Time t (sec.) 
_ Falocity (coef) before (green --> end after cred --) fussy contsol 
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“a 
e 
3 
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° 5 10 15 20 25 0 
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Fig. 6. Plots of the coefficients w11(t), wi1(t) (characteristics of displacement and velocity) before and 
after the control. 


I mm MM 


tN 


i 


He 


! 


= 


lh ii, 


ak arth AH AAA 


Fig. 7. The external and control forces with respect to time. 
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A detailed analysis of estimates for the Fourier coefficients and the rate of convergence of the 
approximate solution are given in the work [15]. The time spectral as well the time collocation 
methods have a high accuracy order of computations, therefore, for small values of N, we can 
obtain good results. Regarding the parameters of the exciting loading forces as many numerical 
experiments have shown changes of the frequency and the amplitude of the harmonic transversal 
loading in the neighborhood of them do not influence much on the control simulation. Therefore, 
ones we compose the FIS for the concrete values of the frequency and amplitude we can use this 
FIS in order to simulate the model for the frequency and amplitude close to the initial ones. 


5. State-space representation for the linearized model 


If we do not consider the nonlinear part in Eq.16 then the system Eq. 16 is a linear model which 
results in a spartial discretization of the Kirchhoff-Love plate model for small deflections 
([11,21,22]). In this case, we can use both “‘S-Function" and *‘State-Space" continuous block 
from the SIMULINK library. 


Let us now write down the linearized model in the state space representation, 


v = Av+ BF, 
w = Cv + DF, (18) 


where v = Ma A= ae Aig B= Pal C= lo HA} D = [0] and F = Pqy + 


Fzy.The size of the matrices A, B, C and D is 2 x N* X 2 x N?. The identity matrix I and 
M, K, C have the size N? x N2. The vectors F and v have the size 2 x N?. 


We obtain analogous representations for the system Eq.17, if we neglect the nonlinear terms. 


A similar algorithm for vibration suppression of the Fourier coefficients for the linear model Eq. 
18 is constructed. In this case, instead of “S-Function”, we employ “‘State-Space" representation 
from the SIMULINK library, which is a widely used tool in the classical control theory. In the 
algorithm for the linear problem in Step 2 we form the matrices for the state-space 
representation. Note, that in the block *‘State-Space", A = A, C = Cand u=F. The initial model- 
diagram for the linearized system is presented in Figure 8 (‘Selector" block corresponds to 
Example 2). The complete model-diagram with control uses the state-space representation inside 
its subsystems. It has a similar structure with the complete model-diagram, Figure 5. The 
“Subsystem 1” and “Subsystem 2” in this diagram are identical with Figure 8. 


sine Wave State-Space [1,10] of 18 


Plots 


Fig. 8. The initial model-diagram for the linear system. 
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Example 2. Here we suppose that the external force Q(t) = sin 25zt is applied at some discrete 
points on the simply supported plate, exactly at the four collocation points: (x,,y,) = (1/(N + 
1),1/(N+1)), (@ry3s) = A/(N + 1),3/(N + 1)), (3,91) = (3/(N + 1),1/( + 1), 
(x3, ¥3) = (3/(N + 1),3/(N + 1)). The physical parameters take the same values as in Example 
1. The number of the basic functions N = 3 and the final time is T = 45. Here q,, = 4,3 = 
931 = 133 = Sin 25nt (see Eq. 15). The initial model-diagram is given in Figure 8 and the 
complete model-diagram has the same structure as in Figure 5. The amplitude in “Sine Wave" 
block is P = (0,0,0,0,0,0,0,0,0,1,0,1,0,0,0,1,0,1)7. The results are presented in Figures 9 and 10. 
The suppression of the first Fourier coefficients is shown. The x-axis (horizontal) denotes t - the 
time of simulation in seconds while the vertical axes denote the first Fourier coefficients in the 
expansions for the displacement and velocity, Eq.3. 
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Fig. 9. Plots of the coefficients w(t), w/1(t) (characteristics of displacement and velocity) before and 
after the control 
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Fig. 10. The external and control forces with respect to time. 


6. Conclusions 


An effective algorithm for vibration suppression of a rectangular elastic plate with an application 
of SIMULINK-MATLAB has been proposed. Two approaches have spatially discretized the 
mechanical plate modelme spectral (Galerkin) method, which was introduced before, and a new 
time spectral-collocation scheme, developed here. The presented algorithm involves model- 
diagrams, created with the use of SIMULINK library and Fuzzy ToolBox. A Mamdani type 
fuzzy inference system has been chosen. 


The techniques have been demonstrated on two examples. We have shown suppression of the 
first Fourier coefficients, which characterize the behavior of the plate. The time spectral 
(Galerkin) method provides more accurate computations. However the collocation method is 
more flexible in the sense of the location of control points. 


Employing of SIMULINK allows quickly and effectively solving the studied control problems in 
linear and nonlinear mechanics and also to improve the presentation of the results. Several 
examples of suppression vibration of a plate with sinusoidal exciting loading forces can be easily 
implemented with the use of the proposed algorithm. The introduced methodology can be 
extended and applied to similar problems. 
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